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Abstract 

The harmonic oscillator as a distinguished dynamical system can be defined not only on 
the Euclidean plane but also on the sphere and on the hyperbolic plane, and more generally 
on any configuration space with constant curvature and with a metric of any signature, either 
Riemannian (definite positive) or Lorentzian (indefinite). In this paper we study the main 
properties of these 'curved' harmonic oscillators simultaneously on any such configuration space, 
using a Cayley- Klein (CK) type approach, with two free parameters Ki , K2 which altogether 
correspond to the possible values for curvature and signature type: the generic Riemannian and 
Lorentzian spaces of constant curvature (sphere S , hyperbolic plane H , AntiDeSitter sphere 
AdS 1+1 and DeSitter sphere dS 1+1 ) appear in this family, with the Euclidean and Minkowski 
spaces as flat limits. 

We solve the equations of motion for the 'curved' harmonic oscillator and obtain explicit 
expressions for the orbits by using three different methods: first by direct integration, second 
by obtaining the general CK version of the Binet's equation and third, as a consequence of 
its superintegrable character. The orbits are conies with centre at the potential origin in any 
CK space, thereby extending this well known Euclidean property to any constant curvature 
configuration space. The final part of the article, that has a more geometric character, presents 
those results of the theory of conies on spaces of constant curvature which are pertinent. 

Keywords: 

The harmonic oscillator. Integrability and Superintegrability. Spaces of constant cur- 
vature. Central potentials. Conies on the sphere and hyperbolic plane. Conies in 
Lorentzian spaces. 

Running title: The harmonic oscillator on spaces with constant curvature. 

PACS codes: 02.30.Hq, 02.40.Ky, 45.20.JJ 

AMS classification: 37J15, 37J35, 70H06, 70H33, 70G65, 

momentum maps, reduction structure of phase space, integration methods 



1 E-mail address: jfc@unizar.es 
E-mail address: mfran@unizar.es 
E-mail address: msn@fta.uva.es 



1 



Harmonic Oscillator on Ricmannian . . . Carincna et al. 



2 



1 Introduction 

In a sense this article can be considered as a sequel or continuation of a previous paper [12] which 
was devoted to the study of mechanical systems on Riemannian configuration spaces with constant 
curvature k ^= 0. Geodesic motion, the theory of symmetries and general results on central potentials 
were discussed in the first part of [12] , while in the second part attention was focused on the Kepler 
problem in S 2 and H 2 . Now, we present a similar analysis for the harmonic oscillator, yet extending 
the scope so as to include also the much less explored cases where the configuration space is a 
Lorentzian manifold with constant curvature. We follow the approach of |12j . which contains the 
fundamental ideas and motivations, and also use the notation, ideas and results discussed in [501151] . 

The study of mechanical systems on Riemannian spaces has been mainly done in connection with 
relativity and gravitation. Nevertheless, before relativity, the study of both Kepler and harmonic 
oscillator potentials in spaces of constant curvature had also been done from the viewpoint of 
classical nonrelativistic mechanics (see [17] for an historical account of the research made until 
the first years of the XX century and references in [12] for more recent papers including also some 
quantum problems on spaces with curvature as, e.g., the hydrogen atom in a spherical or hyperbolic 
geometry). It is interesting to point out that [39], a book on geometry, includes however a final 
chapter devoted to mechanics (the title of this chapter was "Nichteuklidische Mechanik" in the 
original edition but was changed to "Mechanik und spezielle Relativitdtstheorie v in the revised 1923 
edition); in addition to rather general properties linking geometry with mechanics, this chapter 
contains the basics of a study of the harmonic oscillator in constant curvature; polar coordinates 
are used and the approach is basically Newtonian. 

It is well known that the Kepler problem and the harmonic oscillator are the two more important 
supcrintegrablc systems in Euclidean space (see for instance the recent book |16)V and, as it was 
to be expected, they have 'curved versions' which remain superintegrable in spherical or hyperbolic 
configuration spaces |57[ 161] . This known property implicitly underlies some classical papers as 
[56] or [28] . On the other hand, much work has been recently done in the study of superintegrable 
systems in non-Euclidean spaces [H [SI El [HI HDl 1201 HDl [3H HH [351 EH [SD] and this fact has 
intensified the interest for the study of the 'curved' versions of these two systems and their relations 
[3U 1231 EH1 E2] . A further step, which we take in this paper, is to extend these studies also to the 
case where the configuration space is itself a constant curvature Lorentzian manifold. This case 
was definitely not taken into account in the previous papers, and opens some views into a relatively 
unknown field. For some work related to dynamics in Lorentzian manifolds, see [71 l42l [43l 159] . 

The three classical spaces with constant curvature k, to wit, the sphere S 2 with k > 0, Euclidean 
plane E 2 for k — and hyperbolic plane H 2 for k < 0, can be considered as the three different 
instances in the family of homogeneous Ricmannian manifolds V£ = (S 2 , E 2 , H 2 ). A technique 
for considering these three spaces at the same time in a unique family, with the curvature k as a 
parameter k € M was first introduced by Weierstrass and Killing [17] and lies at the origin of the 
so-called Weierstrass model for the hyperbolic plane; for some reason modern presentations usually 
restrict to the standard value k = — 1, thus losing from direct view how some properties depend 
on the curvature. If k is left explicitly, this allows consideration of the 'curved' harmonic oscillator 
(or Kepler problem) on constant curvature spaces as arising from Lagrangians depending on k as 
a parameter and understood as defined in a generic space V* (either (S 2 , E 2 , H 2 )) by using a 
unique K-dependent expression. The 'curved' systems appear as a '^-deformation' of the well known 
Euclidean system and they can be defined without ambiguities because superintegrability picks up 
essentially a unique '^-deformation' among the many (non-superintegrable) potentials having the 
(Euclidean) Kepler or oscillator as its 'k — ► limit'. 
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A convenient tool for the use of k as a parameter are the following ^-trigonometric 'Sine' and 
'Cosine' functions 



as well as the 'Tangent' T K (x) = S K (x)/ C K (x). These functions allow us to write K-dependent 
expressions in a unified way for the whole previous family of spaces V% so that the computations, 
statements and results are unified too. 

This approach, used in some previous papers [3J [T21 H7J EH IE] , has a potentiality beyond the 
unification for the three constant curvature configuration spaces V„ : the mathematically natural 
frame for the 'k as a parameter' idea [TTJ tSH] involves not just one single parameter k, but 

two parameters, Ki and K2, which correspond to a space S K] \ K2 i with constant curvature n\ and 
whose metric is definite (Ricmannian), degenerate or indefinite (Lorentzian) according to the K2 
sign (see e.g. M M 125]). 

In more detail, the plan of this article is as follows: In Sec. II we analyse, in a single run, 
the free geodesic motion on the three Riemannian spaces F 2 = (S 2 ,E 2 ,H 2 ), and in the three 
Lorentzian spaces L 2 K = (AdS* +1 , M 1+1 , dS^ +1 ) as well as the Killing vector fields in the general 
S Kl i K2 \ and its correspondent Noether symmetries. We have divided this section in two subsections, 
using respectively geodesic polar coordinates and geodesic parallel coordinates. 

Sec. Ill, that can be considered as the central part of the article, is devoted to the study of 
the harmonic oscillator on the general space S Kl [ K2 \ ; the specialization to Ki=k, K2 = 1 affords the 
relevant results for the harmonic oscillator in the three 'classical' Riemannian spaces of constant 
curvature (S 2 ,E 2 ,H 2 ), and further specialization to k = leads to the harmonic oscillator in E 2 . 
The results also cover the harmonic oscillator on Lorentzian configuration spaces. 

After an introduction, we have divided this section into three subsections: in the first part we 
consider the equivalent one-dimensional problem and we draw some information and a classifica- 
tion of the orbits previous to obtaining any closed-form solution for the motion; in the second 
part we solve explicitly the problem and we obtain closed expressions of the orbits. The results 
can be obtained by using three different methods: first by direct integration, second by obtaining 
the S ,2 1 [ K2 ] -dependent version of the Binet's equation, and third by exploiting anew the superinte- 
grability of the problem; this last method is closer to the one usually employed in the Euclidean 
oscillator, which leans on its separable character in Cartesian coordinates, a property which is not 
shared by the Kepler potential. The third subsection is devoted to the analysis and classification of 
orbits, with the emphasis restricted to the three 'classical' Riemannian configuration spaces with 
nonzero curvature n. This third subsection also interprets the trajectories obtained above as conies 
in curved spaces; of course all the results reduce to well known Euclidean trajectories — ellipses 
centred in the potential origin — , when we specialize the parameters to their standard Euclidean 
values «i = 0, ft 2 — 1. Our description of conies goes beyond the classical papers on this topics 
published around 1900 (where the emphasis was mainly projective, although some metrical aspects 
are also discussed; see, e.g., [55] and references therein or [T51[3n])- And, at any rate, conies in 
locally Minkowskian were definitely not considered at all in any of these works; probably [5J is the 
first paper dealing with this subject. 

Sec. IV, that has a more geometric character, gives some information on the theory of conies 
in the general CK space S K1 r K2 i , leading to the identification of the harmonic oscillator orbits with 
conies for any value of the parameters m, k.2 (this is, either for Ricmannian as well as for Lorentzian 
configuration spaces), and serves as a geometrical counterpart and complement to the information 




K > 
K = 
K < 



(1) 
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already provided on Sec. III. Additional details are provided mainly in the Riemannian case; we 
plan to discuss the case of a Lorentzian configuration space in more extension in a forthcoming 
paper. Finally, in Sec. V we make some final comments. 

We mention that most results here also hold (in a suitably reformulated way) for the n- 
dimcnsional version of the harmonic oscillator, which is well known to be superintegrable [261 140] . 

2 Dynamics on S^j K y. Geodesic motion, Noether symme- 
tries and constants of motion 

We start by discussing some details on the motion of a particle in a configuration space Sw^i in 
the general CK case, where the parameters «i, K2 may have any real value. This (k\, K2)-dependent 
formalism contains nine essentially different Cayley-Klein spaces, because by scaling the units of 
length by A and of angle by a, the values of K\ and K2 transform as k\ —* A Ki, k,2 —> a 2 K2- 
Then without any loss of generality, any CK space can be brought to its standard form, with either 
K\ = 1,0,-1 and K2 = 1, 0, —I. The standard form of any expression in the CK formalism coincides 
with the result one would obtain by working from the outset in a single space. Should we proceed 
this way, however, the consideration of how details change when there is a variation of curvature 
or signature type would require an additional separate study. The distinctive trait in the CK 
formalism is that the dependence on Ki and n 2 is built-in, and makes a further study of the limiting 
processes fully redundant. For this reason we will keep the general form, with explicit k,\ and «2 
in most of the paper, stressing when required the specific properties holding after a specialization 
for the values of «i, K2- 

The presence of two parameters in the CK family of two-dimensional spaces is related to the 
Cayley-Klein theory of projective metrics, and underlies the length/angle duality which is the 
residue of the general duality in projective geometry when projective metrics are taken into account. 
It is also related to the existence of two commuting involutions in their isomctry Lie algebras [21 j 

Within the generic standard choices K2 = ±1 for the signature type, the CK family S k a R2 \ 
with K2 = 1 includes the three 'classical' Riemannian spaces with constant curvature K\ = K and 
the CK family •S^ 1 [ K2 ] with K2 = — 1 includes the three Lorentzian spaces with constant curvature 
(kincmatically intcrprctablc as homogeneous space-times). In the non-generic case ^2 = the CK 
spaces can be interpreted as the three 1 + 1 non-relativistic space-times, which are limits of the 
spaces with K2 ^ 0; see [23j [44] and references therein. These nine spaces can be conveniently 
displayed in a Table; for more details see [23l [241 [25] . 



Table 1: The nine standard two-dimensional CK spaces S 2 ^ K2 y 





Measure of distance & 


Sign of ki 


Measure of angle 


Elliptic 


Parabolic 


hyperbolic 


& Sign of K2 


Kl = 1 


KI = 


Kl = —1 




Elliptic 


Euclidean 


hyperbolic 


Elliptic K2 = 1 


S 2 


E 2 


H 2 




Co-Euclidean 


Galilean 


Co-Minkowskian 




Oscillating NH 




Expanding NH 


Parabolic K2 = 


ANH 1+1 


G l+l 


NH 1+1 




Co-hyperbolic 


Minkowskian 


Doubly hyperbolic 




Anti-de Sitter 




De Sitter 


hyperbolic K2 = — 1 


AdS 1+1 


M 1+1 


dS 1+1 



On any general two-dimensional Riemannian V 2 or Lorentzian space L 2 , not necessarily of 
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constant curvature, there are two distinguished types of local coordinate systems, 'geodesic parallel' 
and 'geodesic polar', that reduce to the familiar Cartesian and polar coordinates on the Euclidean 
or Minkowskian plane [HI [37] (see Appendix) . The K-dcpendent Kepler problem in V 2 was studied 
in [12] only in polar coordinates but, since the Euclidean oscillator allows separation also in parallel 
coordinates and this property is shared for its 'curved' version, we will use in this paper both types 
of 'geodesic' coordinates: polar (r, </>) and parallel (u,y). 

2.1 Polar coordinates 

The following expression, where ki,k 2 are two real parameters, 

ds 2 = dr 2 + k 2 S 2 Kl (r) d<p 2 , (2) 

represents, in polar coordinates (r, 0), the differential line clement on the space S K r K3 y When 
n 2 = 1, these spaces are, according to Ki >,=, < 0, the three classical Ricmannian spaces V 2 ± = 
, E 2 ,H^ with constant curvature K\. In the three standard cases K\ = 1,0,-1 the metrics 
correspond respectively to the standard sphere S 2 (k\ = 1, K 2 = 1), Euclidean plane E 2 (ki =0, k 2 = 
1), and hyperbolic or Lobachevsky plane H 2 (ki = — 1, k 2 = 1): 

ds 2 | g2 = dr 2 + (sin 2 r) d<f> 2 , ds 2 | E2 = dr 2 + r 2 d<f> 2 , ^s 2 | H2 = dr 2 + (sinh 2 r) def 2 , 

Likewise, when «2 = — 1 these spaces are the pscudo-Ricmannian 2d spaces L 2 X with indefinite 
non-degenerate metric (hence Lorentzian) of constant curvature k% and for the standard values 
Ki = 1, 0, —1 the metrics reduces to: 

ds 2 | Adgl+1 = dr 2 — (sin 2 r) dip 2 , ^ s2 |m 1 + 1 = ^ ~ ^ d4> 2 i ^ s2 |ds 1+1 = ^ ~ ( sm h 2 r ) ^ 2 > 

which correspond to the three standard Lorentzian spaces AdS 1+1 , M 1+1 , dS 1+1 . Here polar coordi- 
nates only cover the region with 'time-like' separation from the origin; unlike the Ricmannian case, 
AdS 1+1 and dS 1+1 are related by a change of sign in the metric, and are thus essentially the same 
space; this transformation is conveyed by the change r *-y ir, which interchanges the regions with 
time-like and space-like separation from the origin. 

The three vector fields Xp x , Jp 2 , Xj, whose coordinate expressions are given by: 

or T K1 (r) d(p 
X P2 = ^S^) - + —. — , 

arc Killing vector fields of the metric, and are well defined for any value of K\,k 2 (even in the 
most degenerate CK space, the Galilean or isotropic space where Ki = 0, k 2 ~0). Each X generates 
a one-parameter group of isometries of the metric and altogether close on a Lie algebra denoted 
«>«i,« a (3): 

[Xj, Xp^\ = — Xp 2 [Xj, Xp 2 ] = k 2 X Pi [Xp x ,Xp 2 \ = — k\Xj. (4) 

Of course, when «i, k 2 are set to any particular values, all these expressions give the pertinent 
ones for the corresponding spaces; this is the trait in all the CK formalism. Notice that only when 
Ki = (Euclidean, Galilean and Minkowskian plane) X Pl and Xp 2 commute. If we restrict to the 
family of classical homogeneous Ricmannian spaces V 2 x with curvature Ki, the corresponding Killing 
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vector fields are given by setting k 2 = 1 in ([3]), and then the angular coordinate appears through the 
circular trigonometric functions in the three spaces, where the radial coordinate appears through 
the Ki-ones, which are either circular, parabolic or hyperbolic according to the sign of k\: 

X Pl \ V 2 = (cos< 



X Pl \ V 2 = (sine 

d_ 

dcf> ■ 



d 


(G Kl (r) . ' 


)- 


dr 


/ d<j> 


d 


fC Kl (r) , 


)~ 


dr 


1 d<f> 



Xj\y2 



Moreover, the Lagrangian for a (free) particle moving in a configuration space S 2 r K2 i is given 
by the kinetic term invariant under the actions of Xp 1 , Xp 2 , Xj arising from the metric: 

C (r,cj), iv,i^) = T^ KuK ^(r,^>,v r ,v^) = 5 (v 2 + k 2 S 2 Kl (r)vl) , 
A general natural Lagrangian in S K j K2 \ (kinetic minus potential term) has the following form 

C(r, <j), v r , V4) = j ( v 2 r + k 2 S 2 Kl (r)v 2 ) - V(r, 4>) , 
in such a way that for k% =0, K 2 = 1 we recover a standard Euclidean system 
L(r, </>, v r , vj,) = lim £(r, 0, w r , w ) = - (v 2 + r 2 vl) - V(r, (j)) , V(r, <t>) = lim V(r, 0) . 

«1— >0,K 2 =1 Z ^ Si- >0,K2=1 

In some particular cases a Lagrangian system can possess the Killing vector fields X Pl , Xp 2 , or Xj 
(or any linear combination of them) as exact Noether symmetries. If we denote by X 1 the natural 
tangent lift to the tangent bundle (velocity phase space) of the vector field X and by 9l the Cartan 
semibasic one-form [5T] 

dL dL 
® L = ~dv~ dr + ~dii~ ^ = Vr dr + K2 ^ ' 

then the basic cases with exact Noether symmetries are the following: 

1. If the potential V(r, <f>) is invariant under Xp 1 , then V(r,<f>) should depend on (r, 0) only 
through an arbitrary function of the single variable z 2 = S Kl (r) S K2 (</>) and then 

P 1 = i(X Pi )6 L = {C K2 {4>)) v r - K 2 (C K1 (r) S«(r) S K2 (0)) V<p 

is a constant of motion. 

2. If the potential V(r,tf>) is invariant under Xp 2 , then V(r, 0) should be an arbitrary function 
of the single variable z\ = S Kl (r) C K2 (</>) only, and then 

P 2 = i(X P2 )6 L = k 2 (S k M) v r + K 2 (C Kl (r) S 2 Kl (r) CM) 

is a constant of motion. 

3. If the potential V(r, 0) is invariant under Xj, then V(r, 0) should be an arbitrary function of 
the single variable r only (V is a central potential) and the constant of motion is: 

J = i(X t J )6 L = K2 S 2 K1 (r)v <i> . 
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Several remarks are pertinent. First, the radial dependence in the momenta P\ 1 P 2 ,J appears 
through fti-trigonometric functions, and hence is sensitive to the curvature; the angular dependence 
is carried through ^-trigonometric functions and in the three classical («2 = 1) Ricmannian spaces 
it appears through cos<^> or sin^, irrespectively of the curvature. Second, the quantities Pi, P 2 , J, 
could be considered as the ordinary linear momenta and angular momentum of a particle moving 
in the configuration space S 2 ^y In terms of these, the kinetic energy P( K i,k 2 ) can be rewritten as 
follows (it is the Casimir of the isometry algebra, [22] ) 

1 n 2 Pl + Pi + K X J 2 



T, 



(«1,K2) 



2 n 2 



showing that, on spaces of (constant) non-zero curvature k±, the angular momentum has a contri- 
bution to the kinetic energy of the system, proportional to the curvature k±. And third, the new 
quantities P 2 and J vanish identically when k 2 = 0. This is linked to the singular character of the 
corresponding Lagrangian, as the metric is degenerate when k 2 = 0. This singular case is however 
not generically singular, but only a very special limit of a regular system, and one may expect 
the geodesic motion to have precisely three non-trivial constants of motion linear in the velocities. 
When working in the general CK scheme, where we want to cover all CK spaces Sw^i (even when 
k 2 = 0) this suggests to consider, instead of Pi,P 2l J, the quantities defined as 

Pi := Pi, V 2 := — , J := — , (5) 

K 2 K 2 

which will be called the CK Noether momenta; when k 2 ^ these are essentially equivalent to 
P\ 1 P 2l J but V\ 1 V 2l J are to be preferred because they remain non vanishing even in the limit 
k 2 — > 0. These will be the momenta used in the rest of the paper; we remark that in the classical 
Ricmannian spaces of constant curvature K\ the three Noether momenta coincide with P\,P 2 , J. 
In 5 , ki[k 2 ] t nc CK Noether momenta are: 

V\ = C K2 (0) v r - k 2 C Kl (r) S Kl (r) S K2 ((j>) v^, , (6) 
V 2 = S K2 (0) v r + C Kl (r) S Kl (r) C K2 (0) , (7) 
J = S^r)^. (8) 

In terms of these Noether momenta, the kinetic energy is well defined for all CK spaces, contains 
always a term V\ and is given by: 

t (kiiK2 ) = \ (Pi 2 + k 2 vI + KlK2 j 2 ) . (9) 

2.2 Parallel coordinates 

The element of arc length ds 2 in the space S 2 < K2 i is given in parallel coordinates (u, y) by 

ds 2 = C 2 KiK2 (y) du 2 + k 2 dy 2 , (10) 
and in the three particular classical Ricmannian standard cases it reduces to 

ds 2 \ s2 = (cos 2 y) du 2 + dy 2 , ds 2 1 E2 = du 2 + dy 2 , ds 2 1 2 = (cosh 2 y) du 2 + dy 2 . 
The Lagrangian of a free particle in S Kl < K2 i has a kinetic term corresponding to the metric: 
C(u,y;v u ,v y ) = T (Kl , K2 )(u,y]v u ,v y ) = \ ( C 2 iK2 (y) v 2 u + ti 2 v 2 y ) . 
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In these coordinates the three Killing vector fields closing the Lie algebra so KliK2 (3) are 

x Pl -- 

Xp 2 - 

Xj -- 

and the associated momenta 

Pi = i{X l Pi )9 L 
P 2 = i(Xp 2 ) Qjj 
J = iX t J )6 L = 

lead to the CK Noether momenta: 

Vi = Gl lK2 (y)v Ul (12) 
V 2 = ki S Kl (u) S KlK2 (y) C KlK2 (y) v u + C K1 (u) v y , (13) 
J = - C Kl (u) S KlK2 (y) C KlK2 (y) v u + S K1 (u) v y . (14) 

Notice how the expressions for the CK momenta specialize in the Euclidean case: 

Vi\- E 2=v u V2\ E 2=v y J\ E 2 = uvy - yv u . 

The general CK expressions can be looked at as a two-parameter deformation of the Galilean 
ones Ki — 0, K 2 — 0, governed by the two constants k\,K2] the Euclidean case is not the natural 
comparison standard in the deformation, as one of the constants is already non- vanishing for E 2 . 
A potential V, now expressed as a function of (u,y), turns out to be invariant under Xp 1 if 

V is an arbitrary function of the single variable y only; this result is simply the translation to 
parallel coordinates of the previous result in polar coordinates — as consequence of the relation 
S KlK2 (y) = S Kl (r) S K2 ((j)) — , so that the auxiliary variable z% turns out to be precisely S KlK2 (y). 
Similar results describe the general form, in coordinates (u, y) of a potential invariant under Xp 2 
or under Xj. 

We close this section with a comment on the expressions obtained for the three vector fields, Xp 1 , 
Xp 21 Xj (for K2 = 1 this point was discussed in |51j). According to the straightcning-out theorem 
PP, a vector field X on a n-manifold V always admits a local coordinate system {x\, . . . ,x n } in 
an appropriate neighbourhood of a regular point X{m) /0,m6 M, such that then it becomes 

Y = ^ k Ck{d/dxk), with Ck = 1, Ck = for k ^ k . We recall that in polar coordinates Xj is 
given by Xj = d/d(j> and now we have obtained that in parallel (u,y) coordinates Xp x takes the 
form Xp 1 = d/du. Similarly, in the complementary 'orthogonal' parallel system (x, v), we obtain 
Xp 2 = d/dv (this parallel system is not used in this paper but is discussed in the Appendix). So, 
these three coordinate systems, (r, <j>), (u, y), and (x, v), are the three appropriated systems (via the 
straightening-out theorem) providing the 'straight' expressions of Xj, Xp l: and Xp 21 respectively. 

3 The harmonic oscillator on S„ \„ i 

The following Lagrangian in the CK space <Sk 1 [ K2 j with curvature n\ and signature type «2 

C{r,4>,v r ^) = \{vl + K 2 $llr)vl)-V HO {r), V HO = T^(r) , (15) 



d_ 

du' 

d d 
= K 1 K 2 S Kl {u)T KlK2 (y) — + C Kl (u) — , 

d d 

= -K 2 C Kl {u) T KlK2 (y) — + S Kl (u) — , (11) 
= Gl lK2 (y) v u , 

= kik 2 S Kl {u) G KlK2 (y) S KlK2 {y) v u + n 2 C Kl (u) v v , , 
-«2 C K1 (u) C K1K2 (y) S K1K2 (y) v u + k 2 S K1 (u) v y , 
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represents the 'harmonic oscillator' in the space S 2 ^^ [501 02]; the potential Vho{t) is 'central' in 
the sense it depends on the radial coordinate only; this dependence involves the label K\ (but not 
k 2 ) and reduces to 



v ho\ k1=1 = \ul tan 2 r , V H o\ K1=0 = V = \ u% r 2 , V H o\ Kl =_i = \ w o tanh 2 r , 

which are the three harmonic oscillator potentials in the three classical standard Riemannian spaces 
S 2 , E 2 , H 2 . The Euclidean function V(r) appears in this formalism as making a separation between 
two different behaviours (see Figure I). In the sphere this potential was considered by Liebmann 
[39] (1905 edition), and later on by Higgs [28] and Leemon [38] . 

Recall r denotes the distance to the origin point computed in the intrinsic metric on the CK 
space. Thus in the classical Riemannian case (V" K 2 ,ft 2 = 1), the potential has a zero (minimum) 
value at the origin and starts growing quadratically with the distance to the origin point, as implied 
by the approximation T K1 (r) ks r around r = 0, which holds for all values of K±. In the flat case 
the potential grows quadratically with any r and approaches an infinite value only when r — > 00. 
When the curvature is non-zero, the behaviour differs in a way depending on the curvature sign. 
When Ki is positive the potential grows faster and tends to infinity at a finite value r = 2 J Ki > 
this is, on the sphere 'equator' (with the origin taken as the pole); the harmonic oscillator on the 
sphere splits the configuration space into two halves by an infinite potential wall on the equator, 
so the spherical harmonic oscillator has two antipodal centres. In the negative curvature case, the 
potential grows slower than in the flat case, and as r — > 00 approaches a plateau, with a (positive) 
finite height Voo := w 2 /(— 2k\). 

The motion in this potential is superintegrable in all CK spaces <S Kl u 2 i since, in addition to the 
angular momentum J , this system is endowed with the following quadratic constants of the motion 

I J 2 = J 2 , 

?2 1 va . /-„ JA 1A) /_ JA _ , ,2 rp2 / \ n 2 



I vl = V\ + W n (r, 4>), Wn (r, 4>) = T 2 Kl (r) C^) , 

I n = VI + W 22 (r, 0) W 22 (r, <f>) = wl T^(r) S 2 2 (<^) , 



I VlV2 = ViV 2 + W 12 (r, <f>) Wia(r, 4>) = ^ T 2 Ki (r) C K2 {<j>) S K2 (<f>) ■ 

Remark that Wn(r, <f>), H ; 22('"j <P)> Wia^ 4>) are we H defined for any CK space, and they do not 
vanish identically in none of them. 

The 'energy' of the motion can be written as: 

I E = \{I n + k 2 I V 2 + K lK2 Ij2) = \{V\ + k 2 V\ + k x k 2 J 2 ) + \u 2 T 2 Kl {r) , (17) 

reducing on E 2 to the known Euclidean expression. In the general CK space S k j K2 ^ the kinetic 
energy is no longer given by (one half of) the 'norm' of the 'momentum vector' V\ + k 2 V 2 but 
contains as an extra contribution the square of the angular momentum, proportional to the curvature 
K\ and thus disappearing in the flat case /ti = 0. 

The four integrals of motion in (|16[) cannot be functionally independent; indeed they satisfy the 
relation: 

l V\ l Vl - (MPs) 2 = ■ 

Taken altogether the constants I-p2 , 1-p2 , 1-p x Vi are the components of a (symmetric) tensor under 
the 'rotation subgroup' SO K2 (2) in any space S K1 u 2 i 



Til Fl2\_ hi Iv 1 V 2 \_ rl + Wn{q u q 2 ) ViV 2 + Wi 2 { qi ,q 2 ) 



T 2X T 22 \IviV 2 Ip* J \ViV 2 +Wi 2 (qi,q 2 ) Vi + W 22 ( qi ,q 2 



(18) 
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Thus the essential property of the Euclidean harmonic oscillator, to have a tensor constant of 
motion (the so called Fradkin tensor [T5J[2H]), survives for the 'curved' harmonic oscillator in any 
CK space S K1 u 2 y This tensor, looked at as the general CK form of the Fradkin tensor, contains 
also a complete set of functionally independent constants of motion, because J 2 = Iji is related 
to the determinant of the Fradkin tensor in a 'universal' way, with a relation explicitly independent 
of Kl, K2 ■ 

det(F) = lu' 2 Ij2 (19) 
Another possible choice for three functionally independent constants is the set {I-p?, I-p?, Ij?} ■ 

3.1 The classification of orbits and the equivalent one-dimensional prob- 
lem 

It is very convenient to introduce in the study of central potentials V(r) on the Euclidean space a 
one-dimensional effective potential V , which governs the radial motion after elimination of the 
ignorable angular coordinate. This procedure allows us to obtain a classification of orbits before 
the integration of the equations [19] . Given a central potential V on the CK spaces S K1 r K2 i one can 
proceed similarly and this leads to the one-dimensional 'effective' potential V cff (r): 



y ett (r) = V(r) + 



2S 2 Kl (r)' 



where the extra term plays the role of the 'centrifugal barrier potential', reducing as it should 
be to J 2 /(2r 2 ) in the standard Euclidean case. Therefore we can classify the orbits and obtain 
some additional information for the harmonic oscillator motion in any CK space S K1 r K2 i simply by 
analyzing the effective potential for the harmonic oscillator motion: 

V cff (r) = \ .1 T» + = \ «g T 2 Ki (r) + -^- } + \^ 2 J 2 , (20) 

which in the three particular one-dimensional problems associated to the standard (/ti = 1, 0, — 1; K% = 
1) sphere S 2 , Euclidean plane E 2 , and Lobachevsky plane H 2 reduces respectively to 

j2 j2 

Vffols* = |^tan 2 r+— ^ = |^ 2 tan 2 r+ ' 2 + \ J 2 , 
a z 2 sin r 2 tan r z 

v&L, - 1^ + ^. 

Vfo| H2 = i^ta n h 2 r + ^ = i. 2 tanhV + ^ ? --ij 2 

The standard Euclidean case k\ = 0, K2 = 1 needs no comment. We discuss in some detail the 
Riemannian cases, allowing any nonzero (positive) values for the curvature Ki ^ and signature 
type K2 > 0. The discussion is made within the generic case J ^ 0, and we will assume J > by 
considering if necessary the reversed motion along the same geometric orbit; for J = the motion 
is actually one-dimensional. From now on we also omit the label HO in the potential when it is 
clear from the context. 

(1) Spherical case: Analysis of the potential V eff for K\ >0, K2>0. 

This corresponds to motion in a sphere, and the natural range for the r coordinate is the interval 
(0, 7r/(- v /Ki), which is the span of r along half a geodesic (half a sphere's large circle). Provided Uq 
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is real as implicitly assumed (the constant in Vho is assumed to be positive) the effective potential 
is always positive, V cff (r) > 0, satisfies the following limits in the boundaries 

lim V off (r) = +00 , lim V cff (r) = +00 , 

r->0 r^7r/(2 v ^T) 

and it has a minimum, whose value we denote E(CF), at the point r m given by T Kl (r TO ) = <Jk%3 /loq 

E{J) := V cff (r m ) = yfewoJ + \k x k 2 3 2 . 

This equivalent potential represents therefore, an asymmetrical well on the 'r line', with two barriers 
of infinite height at r = and r = ir/ (2y / K]"), and one single minimum placed in between (Figure 2). 
Thus, for a fixed value of J the situation is as follows: there is not any possible motion for energies 
E < E{J), there is the (unique) motion r(t) = r m when E = E(J'), and for all energy values 
E > E{^f) the motion of the radial coordinate consists of non-linear one-dimensional oscillations 
between the two 'radial' turning points. On the sphere, the motion with a fixed value of J and the 
minimum compatible energy E = E( J^) is a circular motion, on the circle centred in the origin with 
radius r = r m . The trajectories with E{^J) < E lie in a spherical annulus which always contains 
the circle r = r m \ this annulus grows to cover all the hemisphere when E — > 00. 

All these expressions depend on k,\, k 2 and reduce to the well known expressions for the Euclidean 
oscillator when fti = 0, K2 = 1; notice that the full Euclidean plane appears as the limit K\ — > of a 
single hemisphere, with the sphere's 'equator' going to the Euclidean infinity. 

(2) Hyperbolic case: Analysis of the potential V for K\ <0, k 2 >0. 

This is the hyperbolic Lobachevsky plane case, and here r £ (0, 00). Provided again w>o is real 
as implicitly assumed, the effective potential is always positive, V eff (r) > 0, and has a minimum at 
r = r m whenever T K Jr m ) = ^fi^J j 'loq. It satisfies the following limits in the boundaries: 

lim V off (r) = +00 , lim V off (r) = — 2_ ; 

) — >0 + 00 — 2fti 

thus introducing into the problem a new energy scale, to be denoted Eoo 

Eno : = 



-2ki' 

(remark this scale could be also defined in the sphere case, where its value will anyhow fall outside 
of the physically allowed range of energies [E(J~), 00]). The value of the potential at the minimum is 
given again by E(J) := V cff (r m ), but here the energy scale E^ is placed above E(J): E(J) < E^ 
When K\ <0 the (absolute) values of the hyperbolic type tangent T K .(r) are bounded by , 1 ; 
thus depending on the values of J there are two possible generic situations, according to whether 
J is smaller or larger than an angular momentum scale 

u 



3 o< 



/?C2(— Ki 



1. If J < Joo, this is if ^/k2~J '/ujq < the function V (r) has a real minimum at the 
given by T K1 (r m ) = v / K2 1 //ijJo, with value V (r m ) = ^fn^i^^J + iK 1 n 2 J^ 2 . Motions of the 
radial coordinate with the given value of J will only happen with energies E > E(J). Within 
this range of energies there are two possibilities. For E{J) < E < Eoo the radial coordinate 
is bounded, corresponding to a motion in the hyperbolic plane in an annular region bounded 
by two circles (with radius corresponding to the two turning points in the effective one- 
dimensional radial potential); this annulus reduces to a circle of radius r — r m for E = E(J'). 
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When E — > while E < Eqo the exterior turning point approaches to infinity, reaches it for 
E = Eoo and then the exterior turning point disappears. Therefore for E^ < E the motion of 
the radial coordinate is unbounded, corresponding to motions in the hyperbolic plane which 
stays outside of a geodesic circle centred at the origin, whose radius corresponds to the still 
existing interior turning point. 

2. If J > Joo, this is if ^/k^J jioq> } K the function V cff (r) has not any minimum, and 



decreases monotonically from the potential barrier with infinite height at the origin r = to 
the asymptotic plateau when r — > oo; only energies satisfying E > E^ = u>q/{—2k\) will be 
allowed here and all the trajectories will be unbounded (scattering) open curves. 

These two possible behaviours and its separating case are represented in Figure 3. 

(3) Lorentzian case: Analysis of the potential V eff for k 2 < 0. 

A complete discussion of motion in a Lorentzian configuration space will not be done here. We 
simply mention some traits. First, the 'centrifugal barrier' comes explicitly with the 'opposite' sign 
to the Riemannian one. Second, if r is the distance to the origin, the square r 2 can have any sign 
in a Lorentzian space (with r = on the isotropes through the origin), and thus formally r will 
be real and positive on the region with time-like separation from the origin, but pure imaginary in 
the region with space- like separation. As the harmonic oscillator potential depends on (r), the 
potential will always be real, yet it has the two signs on the two regions with time-like and space- 
like separation from the origin. Hence the harmonic oscillator potential has no isolated minima 
at any proper point, and all points on the isotropes through the origin are extremal, degenerate 
saddle points for the potential function. The possible extremal values for the effective potential 
are still given by the equation T Kl (r rn ) = ^fn^J /u>o. Because of the inexistence of isolated minima 
for the potential, there is no reason to require a positive 0Jq. If we accept a negative sign for 
the constant \uu$, (i.e., purely imaginary values for ivq, the products ujo^/k2 and the quotients 
t/k^/ujo would turn real. In spite of these seemingly strange properties, the Lorentzian harmonic 
oscillator is supcrintcgrablc, and thus completely explicit solutions for its motion can be explicitly 
given, irrespective of any possible physical meanings which will be very different from the familiar 
oscillator. We just remind that the 'inverted' or repulsive Euclidean harmonic oscillator (with 
ujq < 0) is also superintegrable, yet it only has scattering motions along centred hyperbolas, which 
are also conies with a centre at the origin; hence the complete set of conies with centre at the origin 
appear as the complete set of orbits in the harmonic oscillator only if we allow general, unrestricted, 
values, for the strength constants. 

3.2 Determination of the orbits of the harmonic oscillator in S%j K2 i 
3.2.1 Method I: Direct Integration for a central potential 

We have previously obtained, making use of the conservation of the total energy E and of the 
angular momentum J two expressions for r and <fi. Eliminating t between both equations we have 




Jdr 




After the change of variable r — ■+ v with v = 1/ T Kl (r), dr 



S Kl (r) dv, this becomes 
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All this holds for a general V. Particularizing for the harmonic oscillator in S^u^, V{v) = ^lo^/v 2 

3, 2E P cu 2 1 2 

where we have used the notation Ep = E — % k\k 2 J 2 . An integration leads to 



dv 



R(v) 



and the change x — v 2 



'R'(x) 

makes the integration elementary. In this way we arrive at the general solution of the form 

-1—=D-GC K M0-M), (21) 

where of course the constant </>o is trivial in the sense that orbits with different values for <po are 
permuted by isometrics of the space, as corresponds to the central nature of the potential. The 
two constants D, G should be related to the values of the energy and angular momentum E, J as 
it follows from the integration: 

Ep G= J_ 

k 2 J 2 ' n 2 J 



D = T~^2> G^^-^J E 2 -k 2 u 2 J 2 , E p = E-^k 1 k 2 J z . (22) 



Hence the particular orbit corresponding to energy E and angular momentum J is: 



1 E P 



T» n 2 J 2 K 2 J 
which can be also rewritten introducing two new constants A, B as: 

1 C 2 (0-0 O ) , S 2 2 (^-^o) 



W#2> -^U 2 CM-fo)), (23) 



T 2 Ki (r) A 2 B 2 

where 



(24) 



^=D-G -L= K2 (D + G), (25) 
are related to energy and angular momentum by the relations (to be stated later in a simpler form) 

(26) 



1 _ E P - V 'E 2 p - h 2 lo 2 J 2 1 _ E P + V E 2 p - k 2 uj 2 J 2 



A 2 K 2 J 2 ' B 2 K 2 J 2 

Note that in the Riemannian case (when k 2 > 0), then < G < D, and therefore both D — G 
and D + G are positive, thus making well adapted the notation we have chosen for A 2 ,B 2 ; in the 
Lorentzian case things behave differently. 

All these equations apply for any CK space. In particular, in the three classical standard 
Riemannian spaces (S 2 , E 2 , H 2 ) with curvature k\ = 1, 0, —1; k 2 = 1, and chosing the origin for so 
that 4>o = 0, the orbit equations are: 

— 1 — — — cos 2 6 + -L sin 2 6 in S 2 , 



tan 2 r A 2 B 2 



— — — cos 2 6 + — r sin 2 6 in E 2 



A 2 B 2 



1 1 o / 1 • 2 i tt9 

- cos m A sin m H 



tanh 2 r A 2 B 2 
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Let us close this direct integration approach with three observations. First, it turns out that 
in any of the nine CK spaces, the HO orbit is always a conic with centre at the potential origin, 
where 'conic' has to be taken in a metric sense, relative to the intrinsic metric in each space. 
This follows from the geometrical study to be presented in Section 4. In the Euclidean case the 
quantities A and B are directly the ellipse semiaxes; in the case K\ 7^ the semiaxes (understood 
in terms of the intrinsic metric) are related to A, B by some relations which will involve ki,k 2 . 
Second, the connection between the coefficients D, G or A, B and the energy and momentum, when 
reexpressed in terms of the quantity Ep does not depend on k\. Hence in the sphere or in the 
hyperbolic plane this relation has the same form as in Euclidean space; something similar happens 
in the Kepler problem in Ricmannian curved spaces |12j . Third, both the method and the results 
obtained show a close similarity with the Euclidean ones. The classical and well known change of 
variable r^v = 1/r admits as a generalization the K^-dcpcndcnt change r—*v = l/T K1 (r) which 
affords a significant simplification for all values of k%; of course this change reduces to the Euclidean 
one r — > v — 1/r for Ki = 0. 



3.2.2 Method II: Equation of Binet 

The expression (|5J) of the angular momentum J determines a relation between the differentials of 
the time t and the angle <f> 

Jdt = S 2 Kl (r)dct>. 

The corresponding relation between the derivatives with respect to t and 4> is 

d ( J \ d 



so that the second derivative with respect to t is given by 



dt 2 



J 



J 



Introducing this in the radial equation r = n- 2 S K1 (r) G K1 (r)4> 2 — 4^, it becomes 



( 



J \ d J-/ J 

S 2 Kl (r)Jd4vs 2 Ki (r) 



J \dr 



J 2 



dV 

dr 



This equation can be simplified in two steps: the term in brackets in the l.h.s. can be rewritten by 
making use of 

( J \ dr _ d ( 1 



to obtain 



-J' 



d 2 , 



-j- 



= K 2 J 2 



T Kl (r), 
1 dV 



b 2 \T K1 (r)J T B1 (r) dr ' 

and then we introduce the change r — > v with the potential V(v) considered as a function of v. In 
this way we arrive at the differential equation of the orbit 



12 + K2V = 



1 dV 



that permits us to obtain <j> as a function of v for the given potential considered as function of v: 

-(1/2) 



fc = ±/{c-(^)v 



dv . 
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(A ± sign is pertinent if both signs of the angular momentum along a given orbit are considered; 
notice this equation coincides with the one obtained in the previous subsection) Let us now par- 
ticularize for the harmonic oscillator V = ^lu'q (l/v 2 ) in the CK space S^^y. the equation itself 
reduces to a nonlinear equation of Pinney-Ermakov type 



d 2 v Lu 2 1 



whose general solution, further to the parameter J already present in the equation depends on two 
independent integration constants D, cf>o and has the form 



v = ^D-GC K M4-M), G = JD 2 °— (27) 

which coincides with the general orbit obtained before. The differential equation of the orbit for 
the variable v , usually known in the Euclidean case as Binet's Equation, is essentially preserved by 
the Ki, ^-deformation. Indeed, for the three classical Ricmannian spaces (k 2 = 1) the new variable 
v = 1/r deforms to v = l/T Kl (r) but the equation by itself remains invariant. 

3.2.3 Method III: Superintegrability in parallel coordinates 

In term of the parallel coordinates (u, y) the Lagrangian which represents the motion of a particle 
under an harmonic oscillator potential V in the CK space S k j K2 \ is: 

C(u, y, v u ,v y ) = | {G 2 KiK iy)v u 2 + k 2 v 2 ) - \ w jj V(u, y) , V(«, y) = + n 2 T 2 KiK2 (y) (28) 

The expression of the potential, which correspond to the function on S K] \ K2 i given previously in 
polar coordinates, displays its separability in (u,y) coordinates. The potential reduces to 

/ tan 2 u 2 \ i 2 ( tanh 2 u 2 

Y~ + tan y) V H 2 = j 72 + tanh ! 

cos 2 y J " 2 u V cosh y 



in the two particular cases of the standard unit sphere S 2 and Lobachcvsky plane H 2 , and to 

V| E2 = V = | u 2 (u 2 + y 2 ) = \ cj 2 (x 2 + y 2 ) 

in the Euclidean case (where we recall the equality u = x for n\ = 0). 

Since we have already solved the dynamics in polar coordinates (r, <p) > we can make use of the 
expressions relating parallel with polar coordinates (where the positive 'it axis' y = is taken to 
coincide with the polar axis = 0; see Appendix). Then the orbit equation (taking 0o = to avoid 
inessential complications) 

(e p ^E 2 p -k^U 2 ) Tl(r) C 2 K2 (cj>) + (e p + ^ E 2 P - n 2 u 2 J 2 ) T 2 x (r) C^ 2 (0) = k 2 J 2 , 
becomes when written in coordinates (u, y) 

(Ep ^E 2 p -n 2 u: 2 J 2 ) T 2 » + (Ep + ^El-^J 2 ) ( ^^)" = ^J 2 , 

or 
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The coefficients A, B are related to the values of the constants of motion. By using the identity 

- + k 2 T KlK2 (y) = T Ki (u) + K 2 — 2 — 



the energy constant 

<(«) 



can be rewritten as a linear combination of the constants related to the supcrintcgrability of the 
harmonic oscillator: 

I E = -(Jpa + K 2 /p| + K 1 K 2 I J 2) 



where the quadratic constants of the motion specific to the harmonic oscillator in the space S 2 _^ 



are: 

^ = ^ 2 



= P 2 2 + W22 (r, <t>) W 22 (u, y) = cu 2 ^f^- ( 31 ) 

As in polar coordinates, all the functions Wn(u, y), W 2 2{u, y), Wi2(w, y) are defined for any CK 
space, and they do not vanish identically in none of them. In the Euclidean case, the CK momenta 
Pi,V 2 and the constants of motion reduce as they should to: 



\\ E 2=V U T 2 \ E 2=V y I V 2 



vI+ujIu 2 , / Pi p 2 | E 2= v u v.y + Luluy, I V 2 



2 , 22 
v y + uj y . 



E 2 

(32) 

where once more we recall that for K\ =0, we have u = x. These constants are the elements of the 
Fradkin tensor [181129] . 

Back to the CK general case, and defining: 



2E 1 :=I V 2=Vf+u J 2 Tl i (u) (33) 

T K!K 2 (y) 



2E 2 ^I V 2=Vl+u J l^j4- ( 34 ) 



the energy can be written as 

E = E 1 + k 2 E 2 + \h x k 2 J 2 

where E\,E<i are the curved analogues to the 'partial' energies associated to the one-dimensional 
harmonic motions whose linear superposition provides the most general 2d Euclidean harmonic 
oscillator motion; remind that a curved configuration space is not an affine space, thus strictly 
speaking there is no a well defined way to superpose motions. The value of the remaining constant 
of motion I-p\Vi is equal to zero for orbits with 4>q = (this choice amounts to diagonalize the 
Fradkin tensor, and 2E\ 1 2E2 are the eigenvalues). In this case, the relation (|19[) gives 

AE X E 2 = oj 2 J 2 , 

henceforth implying 

{Ei + K2E2) 2 - {-Ei + k 2 E 2 ) 2 = Ak 2 EiE 2 = k 2 ujIJ 2 . (35) 
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The quantity Ep introduced before reduces simply to the 'sum' of the two E\, E2 contributions 

Ep = Ei + K2E2 

and some simple algebra leads to neater expressions for the relations among physical constants and 
conic coefficients: 

E 1 + k 2 E 2 -E x + k 2 E 2 2 2Et u2 2E 2 

U = — — , G = -2 A =—5- B =—3-. (3b) 

K 2 J Z K 2 J Z LU$ W5 

Direct observation in the expressions for the 'supcrintcgrability constants' leads also to the 
following identity, holding for all S 2 ^ K .^ : 

T 22 T 2 Ki (u) - 2T 12 T K1 (u) ( ) + Fu ( ) ' = n 2 J 2 . (37) 

This property is proved by direct computation, and can be interpreted as a relation which must 
hold between the coordinates («, y) along a given motion, so this is precisely the orbit equation, 
which can thus be obtained directly from the superintegrable character. This is the most general 
orbit, but if the coordinate axes are chosen so that the Fradkin tensor is diagonal, then = 

and the equation coincides, of course, with (|29[) . In the particular Ki=0, k 2 = 1 Euclidean case we 
obtain 

T 22 x 2 - 2Ti 2 xy + Tuy 2 = J 2 
that, as it is well known, represents an ellipse in the Euclidean plane E 2 . 

3.3 The period of the harmonic oscillator 

The Euclidean oscillator is the classical example of an isochronous system, with the same period 
T = 2tt/loq for all orbits. A natural question is whether or not the 'curved' oscillator also inherits 
this property. We discuss this problem for the closed orbits in Riemannian configuration spaces, 
where n 2 > 0. The strategy consists in choosing the orbits with 0o = 0, and computing the time 
spent by the particle between the points = and <j) = 2 J Ki> , which will always be reached along 
a closed orbit. By symmetry reasons this time equals one fourth of the orbit period. 
The angular velocity, taken from the angular momentum first integral: 

; d(j> J 



dt S^r) 

can be rewritten by successively using the identity 

1 1 



S^(r) T 2 Kl (r) 



and the orbit equation 



^ = ^C 2 2 (0) + ^S 2 2 (0), 

as 

JM^2E- 1 C ^ )+ 2E 2 S ^ <t>)+h ^ 
where by introducing double angles and simplifying we get: 

An 2El E 2 d^ M = / + + + / _ Ei + K2E2 



Juj 1 dt dt 
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By integrating along one fourth of a complete closed orbit: 



T = 4 / dt = 4n 2 J 
Jo Jo 



a + P cos(2y / K20) ' 



where a and (3 denote the following expressions: 

a := (Ei + k 2 E 2 + k iK2 J 2 ^ (3:=(-E t + k 2 E 2 
The change £ = 2 V / K20 makes the integration elementary 

T = 2JTT 2 J / = 2JJ^J- 



Jo a + /3cosC v ^Jo? - 1 
Now some algebra and the use of the relations (|35|) leads to 

2tt 1 
T= - 



1k x E 



an exact result whose expansion in powers of the curvature has the correct Euclidean period as the 
zero-th order term, with a fractional first-order correction by the dimcnsionlcss quotient k\E /uj^: 



2tt/ E \ . 2tt 

T=—(l + ), T\„ = — 



Thus when the configuration space has curvature the period ceases to be the same for all orbits 
(as it was to be expected), yet it depends only on the total energy; this property is characteristic for 
all closed orbits of superintegrable systems. On the sphere, the period tends to infinity only when 
the orbit approaches the equator, where E — > oo, but on the hyperbolic plane one might expect the 
period to diverge precisely when the orbit character changes from a closed orbit to scattering open 
orbit. From the previous section we know this happens for E = E^ and this is precisely the value 
which makes 1 + 2k \ e = 0, and hence T = oo. For energies above E^ the derivation as provided 
is not strictly applicable (because in these cases the orbit does not reach any point with (f> = 2 ^J K2 ; 
this is reflected in the formally imaginary result got for T in these cases). 



4 Conies in spaces of constant curvature 

In this section we give a brief geometric description of conies in the CK spaces, focusing on the three 
classical Riemannian (k 2 > 0) constant curvature spaces S^ 1 ,E 2 ,H^ 1 and emphasizing only those 
aspects relevant in relation with oscillator motion in these spaces. This description is intended to 
be self-contained, but as far as conies in these spaces are concerned, it is also complementary to 
comments made in our previous paper devoted to Kepler problem [12] . 

We start by recalling the basics. Any S k t K2 i is a homogeneous symmetric space, and any 
homogeneous symmetric space has a canonical connection which is always unique and well defined. 
Restricting to the homogeneous spaces S K r K2 i, and due to the quasi-simplicity of the involved Lie 
algebras so KliK2 (3), there is always an invariant metric g in S^j^i which is non-degenerate whenever 
k 2 j^0. The canonical connection is always compatible with the metric. Generically, i.e., for n 2 y^0, 
we have more: the canonical connection of S 2 ^ K2 ^ as symmetric space coincides with the Levi- 
Civita (metric) connection associated to the metric. In the degenerate case k 2 = however, there 
are many connections compatible with the (degenerate) metric, and the canonical connection of 
the homogeneous symmetric space is singled out among the many connections compatible with the 
metric. 
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By lines will mean the autoparallels of this canonical connection. Generically (when «2 7^ 0), 
these coincide with the extremal curves of the length functional. In the Lorentzian case there are 
two generic types of geodesies: time-like, with real length and space-like, with 'imaginary' length 
relative to the main metric g (alternatively one may think of another space-like length along space- 
like curves, corresponding to the 'companion metric' g/«2; the space-like curves have a real 'length' 
as computed relative to this metric yet this choice will not be used in this paper). Thus from now 
on, all distances between points or between a point and a line will refer to lengths along geodesies 
measured relatively to the main metric (hence when K2 < distances might be either real, vanishing 
or pure imaginary). 

The geometric definition of conies, which makes sense in any 2-d space of constant curvature K\ 
and non-degenerate metric (/t2 7^0) involves focal elements, i.e., either points which arc assumed 
oriented, or lines which can be oriented and co-oriented. In any such space, and by definition: 

An ellipse/hyperbola will be the set of points with a constant sum/difference 2a of distances 
ri,r2 to two fixed points F\,F2 separated a distance 2/ and called foci. 

A parabola will be the set of points with a constant sum/difference of distances n, r-i to a fixed 
point F\, called focus, and to a fixed line fi, called focal line, (f2 is assumed to be oriented); the 
oriented distance between F\ and fi plays here the role of focal separation. 

An ultraellipse/ultrahyperbola will be the set of points with a constant sum/difference 2a of 
oriented distances f\ , f% to two fixed intersecting lines f\ , fi separated by an angle 2F and called 
focal lines. 

In the generic CK 2-d space of constant curvature K\ ^ 0, K2 ^ these three pairs of curves, each 
pair sharing the same focal elements, are the generic conies, and all the remaining possible conies 
are either particular instances with focal separation vanishing (/ = 0, cf> = 0, F = 0) or limiting 
cases, where some conic elements go to infinity (if possible at all); both particular and limiting 
cases can be obtained as suitable limits from the generic conies. 

An important observation is that the Euclidean plane is not generic neither among the complete 
family of CK spaces, nor among the restricted family of Ricmannian constant curvature planes. 
Thus some Euclidean properties of conies are very special and do not provide a good viewpoint 
to assess the ki 7^ properties. For instance, in the hyperbolic plane a given ultracllipsc, 
determined according to its definition by a pair of intersecting focal lines, has not only this pair 
of focal lines, but altogether three such pairs, the remaining two pairs being the non-intersecting 
pairs of focal lines obtained by joining in the two possible ways the endpoints of the (intersecting) 
initial pair. The distances between the two members of each pair plays the role of focal separations 
in such cases; and are in a suitable sense complementary (precisely, the paralelism angles of half 
these distances are complementary in the ordinary angular sense). A further important detail is 
that in , ultrahypcrbolas for a given pair of non-intersecting focal lines are ultraellipses for the 
matching pair of non-intersecting focal lines (with the same end-points), hence any ultraellipse can 
also be understood as an ultrahyperbola, unlike ellipses and hyperbolas in . 

To make contact with the results in the physical part, we now draw our attention precisely to 
the conies with centre, because the orbit of a particle in the oscillator potential has a centre 
of symmetry at the origin of the potential. In any CK space ellipses/hyperbolas have centre, as 
do ultracllipscs/ultrahypcrbolas as well as their common limits when the 'major' scmiaxis go to 
infinity (if possible at all; this might happen in the hyperbolic plane, where these limiting conies are 
cquidistants). Further analysis shows that conies appearing in the oscillator problem are generically 
ellipses (but not hyperbolas) and ultraellipses, with circles as the 'equilateral' case of an ellipse 
with equal semiaxes, and eventually cquidistants as the non-generic limiting orbits. Parabolas have 
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however no centre. This mean that oscillator orbits cannot be parabolas in any CK space, and 
henceforth these conies will be disregarded from now on. 

In the sphere S 2 X there are no either points nor lines at infinity, thus there are no limiting 
cases in the sense they will appear in . This means that on the sphere there is a single type of 
harmonic oscillator orbits: gcncrically ellipses on the sphere, with two limiting cases, lines through 
the origin (for vanishing angular momentum) and circles centred in the origin. 

In the Euclidean plane E 2 ultraellipses with any focal angle F are just pairs of parallel lines, in 
directions parallel to the two bisectors of the focal angle; in the standard description of Euclidean 
conies these are degenerate, and as far as harmonic oscillator trajectories these are unphysical and 
correspond to infinite energy and angular momentum. 

The hyperbolic case H 2 . with negative constant curvature is more interesting. The limiting 
conies relevant to the oscillator problem are obtained when the major axis of the ellipse tend to 
infinity, while the minor axis remains finite; the ellipse tends to a limiting conic, which turns out to 
be a equidistant curve to the major axis; this equidistant is also the limiting conic obtained from 
an ultraellipse when its 'major axis' tend to infinity. In E 2 these limiting conies collapse precisely 
to pairs of parallel lines, but here they appear for finite values for E and J . 

4.1 Analysis of the orbits 

4.1.1 The general CK configuration space 

In the physical part (Section 3) we obtained, in any S^r^i, the equation of the orbit (|24p. which 
contains three free parameters. One of them, <f>Q, can be set to zero by taking appropriately the 
direction for the origin of angles. Thus the orbit equation is: 



1 ClU) , S 2 2 (</>) 



T 2 Ki (r) ^ 2 5 2 ' 

which depends on two constants A and B. In the Euclidean case, this equation reduces to: 

1 _ cos 2 (0) sin 2 ((f)) 

^2 A 2 + B 2 > 



(38) 



(39) 



which is an ellipse with semiaxes A,B. Again in the general CK space, the 'physical' constants 
E, Ei, Ei, J are related with the parameters A, B by the relations 

A 2 = ^, B^ 2 -^, A*B> = £. (40) 

UJ Q UJ Q U! Q 

The type of the orbit as a conic in 5^i[k 2 ] depends on the space itself (the values of Ki,k 2 ) an( l 
on the constants E, E\, E2,J, but the explicit form of the relations among E,E\,E 2 ,J and A, B 
turns out to be independent of the CK parameters k\,k 2 ; for the Euclidean case K\ =0, k 2 = 1 they 
reduce to the well known expressions, and in this case the coefficients A, B in ([39)1 are directly the 
orbit semiaxes. The total energy and 'partial' translational energy are given by 

E = E 1 +k 2 E 2 + ±k 1 k 2 J 2 , E p = E 1 + k 2 E 2 . (41) 



The last physical constant is not the total energy, though it coincides with the total energy in all 
flat configuration space); it can be considered as a kind of 'translational' part of the energy. 
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4.1.2 Orbits in a Riemannian CK configuration space 

We now discuss in more detail the standard Riemannian case (k 2 = 1), with special emphasis in the 
two nonzero curvature cases. The type of the orbit depends on the values of the constants A, B. We 
consider altogether the family of all conies whose equation, using polar coordinates in S 2 ,E 2 ,H 2 , 
with the focal symmetry axis of the conic orbit as <f> = 0, may be written in the form: 

1 _ cos 2 Q) sin 2 (0) 
T 2 i(r) A 2 + B 2 • I J 

with A > B > 0. When K\ > this equation describes spherical ellipses, with spherical circles 
as the particular case A = B and with a spherical line (the equator) when A = B = oo. When 
fti < 0, and depending on the values of A, B, this curve which is always a conic might be generically 
cither an ellipse or an ultracllipse in H 2 . . Before discussing each case separately, notice that in 
the three Ricmannian cases the dependence on the polar angle <f> is exactly the same as in the 
Euclidean case. The quantities A, B are related to the lengths of the two semiaxes, but as one 
could expect, the details depend on the sign of k%: if Ki > the range of values of T rei (r) is the 
whole real line (completed with oo), but when n\ < 0, the values of T rei (r) are confined to the 
interval [—l/y/—Ki, l/y/—Ki\, which reduces to [0, l/y/—Ki\ for positive values for r. The values 
for r at <fi = and <p = n/2 will play the role of major a and minor b semiaxis of the conic; then 
the relation between the constants A, B and the geometric semiaxes is 

T K1 (a) = A, T K1K2 (6) = B . (43) 

These relations identify the geometric meaning of the constants A, B. In the non-negative curvature 
case, «i > 0, any A, B will determine uniquely a, b, because the range of the circular and parabolic 
tangent is the whole real line. But this is not so when K\ < 0, for then this equation will define 
real values of a, b only when A,B< 1/ ^J—k\. If both A, B are larger than l/y'— ki, then the curve 
determined by the equation (|42[) is empty (as the values in the r.h.s are never got for any real 
value of r). Hence, we shall always assume B < l/^/—m, so the only alternatives to study in the 
hyperbolic plane H 2 , according to A is smaller than, equal to or larger than 1 / \f—K\. 

We remark that in the Riemannian case the origin of angles can always be chosen so that A > B, 
and the major axis is precisely the focal axis. We shall assume this choice; remark however that in 
the Lorentzian case < there would be two cases to be discussed, in agreement with the existence 
of two kinds of separation. 

• (Standard) Spherical space k\ >0; K2 = l. The polar equation of the orbit is: 

K\ cos 2 <p sin 2 (f> cos 2 <p sin 2 <f) 

tan 2 (^r) =_ 4^ + ~B^ = lt>) + Tl>)- {U) 

For any values for A, B with A > B, the values a, b (with a > b) are uniquely determined, and 
belong to the interval [0, 7r/(2y / K]")]. This curve is always closed, and is a spherical ellipse with 
centre at the potential origin. If this point is taken as the 'pole', the complete orbit is contained in 
one of the half-spheres bounded by the 'equator', because the r.h.s. of eq (|4"4"|) is bounded by below, 
so the value r = 7r/(2 v / 7cT) (for which x^.V) = ^ cann °t be attained. For a fixed value of J the 
minimal value of the total energy corresponds to circles whose radius r m satisfies T Kl (r m ) = J~/loq 
and have total energy E cir = E{J) := u J + ^KxJ 2 . The values of the energies of the possible 
motions for a given J lie in the interval [E(J'), oo]; with our choice for origin of angles Ei > E 2 - 
Any oscillator orbit through any point in the upper/lower half-sphere centred at the origin is always 
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completely contained in that half-sphere, and from that viewpoint, there is a single limiting orbit, 
the 'equator' which corresponds to r m = it/(2^/~ki) and a = b = 7t/(2 v / k7) and A = B = oo; the 
physical constants for this orbit are J = oo, E = oo. 

• (Standard) Hyperbolic space K\ < 0; K2 = 1- Let us now consider the general case of arbitrary 
negative curvature. The orbit equation is: 

—K\ cos 2 (j) sin 2 

tanh 2 (7^r) = ~A^~ + "b 2 "' ( 5) 

This curve is always a conic in the hyperbolic plane, but its type depends on the values of A, B. 
This was to be expected from the discussion for the effective potential, but it is worthy to look at 
this situation from a more geometrical perspective. 

Whitin the family of conies (|4"5)) one can expect ellipses with semiaxes a, b. These quantities 
should be related with A, B by means of ([43]) . This is only possible when A 2 < 1/ — n\, B 2 < 1/ — «i, 
for only if this condition holds the values A, B can be written as hyperbolic tangents of actual values 
a, b. ranging in the interval [0, oo]: 

T 2 K1 (a)=^ 2 , T 2 K1 (b) = B 2 . 
Its equation in polar coordinates is: 



-K\ cos 2 (ft sin 2 1 



tanh 2 (7^r) T^(a) T 2 Ki (b) 



(46) 



What happens when A 2 or B 2 lie out of the former range? If both A 2 > l/(— B 2 > l/(— 
then for any (j) it would follow tanh 2 ( v / — Kir) > 1 a condition which cannot be satisfied in any 
proper point of H 2 . and this situation cannot produce any oscillator orbit. Then we are left with 
the case A 2 > l/(— «i), -B 2 < l/(— In this case the conic (|45j) is not an hyperbolic ellipse, 
but a ultraellipse, and formally its semiaxis is not an actual distance. To cater for these cases we 
introduce another real quantity a, formally complementary to the would-be semiaxis a, and which 
is related to A by 

1 f =A. 

—K% T Kl (a) 

Notice that when a £ [0, oo], the function l/(/cf T re (a)) ranges in the interval [1/ — «i,oo] and 
therefore, when considered altogether with T 2 i (a), the union of the ranges of the two functions 
(a) and (o)) fills the real line. In the following we will refer to a as the ultraellipse 

'semiaxis'; the polar equation of the ultraellipse with semiaxes a and b is: 

_n, = jos 2 ^ ^n^ = ^n 2 ^ 

tanh 2 (^ r) T 2 1K2 (6) 1 ^ ' T 2 iK2 (6) 

This type of conic in H 2 . has no generic analogue in the Euclidean plane; the Euclidean 'limit' 
K\ = of (|47p is a straight line, which would correspond to an Euclidean oscillator orbit with infinite 
energy and angular momentum. On the hyperbolic plane these oscillator orbits reach the spatial 
infinity with finite values for energy and angular momentum. 

Therefore, in H 2 , the equation of the complete family of conies we are considering is given, 
generically, by one of the two mutually exclusive possibilities: 

-Ki cos 2 (j) , sin 2 -ah 2rp2/~x 2 , , sin2 ^ U q\ 



tanh 2 (^Tr) T^a) T^(6) tan^(V=^r) T^(6) 
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where b G [0 < b < oo] and either 0<a<°ooroo>a>0. Both types of orbits have a common 
limiting case, either a — > oo or a — * oo, which corresponds to the conic 

— Ki cos 2 d> sin 2 (A , s9 , sin 2 d> 

tanh 2 (v^r) ^ T 2 iK2 (6) W T 2 ^) 

This conic is neither an ellipse nor an ultracllipsc in the hyperbolic plane, but can be obtained from 
either of these types as the limit a — ► oo or a — » oo while b is fixed. In the Euclidean case, the 
limiting ellipse is clearly a straight line, equidistant to the major axis and with b as equidistance. In 
the hyperbolic plane, this conic is an equidistant curve, with equidistance b from the a;-axis. This 
means that in the case of the hyperbolic plane H 2 the family of conies we are considering includes 
conies intersecting the major axis at proper points (the ellipse), at only at infinity (the equidistant) 
or not intersecting it at all (the ultraellipse). 

Of course the two expressions (|48p can be used as well when n\ > but then each of the two 
alternatives covers actually all cases. And further, only the first possibility has a sensible Euclidean 
limit because the k± — ► limit of the a family lies completely at the infinity; the existence of this 
family as a different set of conies is specific to the hyperbolic plane. 

In Section 3 we discussed the appearance of some energy and angular momentum standards 
^oojiToo m the hyperbolic plane. As it might be expected, these correspond to the transition 
among different types of oscillator orbits. To be precise, very simple calculations show that in H 2 : 

• All elliptic orbits have energies E < Eoc and angular momentum J < 

• The equidistant orbits have energy E = E^ and angular momentum J < CToo 

• The ultraclliptic orbits have energies E > E^ and the angular momentum J can have any value. 
The ultraellipses with J < are those with T Kl (6) < T Kl (a), while those with T Kl (6) > T Kl (a) 
have J > Joo. Separating these behaviours are the 'equilateral' ultraellipses with T Kl (6) = T Kl (a), 
all of which have J = Joo . 

• The 'largest circular orbit', with radius r = oo is a common limit from elliptic orbits when 
a — ► oo,6 — > oo or from equidistant orbits, when b — > oo. This orbit has energy E = E^ and 
angular momentum J = 

It is clear that when the total energy E is smaller than Eoo, the trajectory (an ellipse) is bounded 
and the motion is periodic, while for E equal to or larger than Eoo, the motion is not periodic and 
the orbit goes to spatial infinity. 

We finally give the relation between the two 'geometric' quantities, the ellipse (or ultraellipse) 
semiaxes a (or a) and b, the parameters A, B appearing in the canonical form of the orbit (|38p and 
the angular momentum J and energies E, E\,Ei\ 

^T 2 Ki (a)=^A 2 = 2E 1 , or cj 2 1 =wlA 2 = 2E u (50) 

ulTl iK2 {b)=u 2 B 2 ^2E 2 , (51) 

J = cj T Kl (a) T KlK2 (b), or J = uj * T KlK2 (b), (52) 

-ki I Kl (a) 



i. .2 



or E = 1 ^{ W ^^+^^)>. (53) 
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5 Final comments and outlook 

We have discussed and completely solved the harmonic oscillator problem simultaneously on the nine 
2d spaces of constant curvature and metric of either signature type. As stated in the introduction, 
one of the fundamental characteristics of this approach is the use of two parameters n\, «2 in such a 
way that all properties we have obtained reduce to the appropriate property for the system on each 
particular space. Generically these spaces include either the sphere S 2 , the hyperbolic plane H 2 , 
the De Sitter dS 1+1 or the AntiDeSitter spaces AdS 1+1 when the corresponding values of K\, K2 are 
appropriately set. Euclidean or Minkowskian spaces arise as the very particular (but important) 
flat cases k,\ =0 with a non-degenerate metric K2 7^0. So, we can sum up the results pointing out 
some important facts: 

• The harmonic oscillator is not a specific or special system living only on the Euclidean space. 
In any of the nine CK spaces S%^ K ^ there is a system with all the outstanding properties of 
the Euclidean harmonic oscillator, worth of the name 'curved harmonic oscillator'. 

• Motion in the 'curved harmonic oscillator' potential on any CK space is superintcgrable. 

The orbits of motion in this potential are always conies with centre at the origin of potential, 
in all the CK spaces and in the sense of the intrinsic geometry in S^ K2 y These conies always 
include ellipses, with straight lines trough the origin as limiting cases for J = 0, but other types of 
conies might also appear in some spaces. For the three Riemannian spaces of constant curvature, 
the situation is as follows: 

• In the sphere case, there is a single type of oscillator orbits; all such orbits are ellipses and 
are confined to an hemisphere centred in the oscillator centre; a 'limiting' case is the 'equator' 
orbit, which appears as an ellipse with major axis equal to half the equator's length (on the 
sphere this condition implies both axis to be equal); this orbit has infinite energy and angular 
momentum. 

• In the Euclidean case (the n\ — > limit of the sphere case), all oscillator orbits are ellipses, 
and a limiting case happens when the orbit tends to a straight line not through the origin, 
with infinite angular momentum and energy. 

• In the hyperbolic plane, orbits with values of E smaller than E x are ellipses, and for them 
J should be smaller than Jqq. This ellipse reaches the spatial infinity when E = E^; this 
happens when the focal distance goes to infinity but the minor axis stays finite, and the 
limiting curve is an equidistant, with base curve the major axis and equidistance distance 
equal to the minor semiaxis. For values of E > E^ the orbit is (a branch of) an ultracllipsc, 
a conic reaching the spatial infinity, and in this case all values of J (even larger than j7oo) are 
allowed; in the limit where E, J go to infinity, the orbit tends to a hyperbolic straight line, 
which appears here as a particular limiting case of an ultracllipsc. 

Then, as compared with the Euclidean case, the new trait appearing in the hyperbolic plane 
is the 'splitting' of the Euclidean 'last' singular straight line in a full family of orbits: a 'first' 
orbit reaching the spatial infinity, then a full family of ultraelliptic orbits and finally a straight line 
orbit (Figure 4). And the new trait in the sphere case is that the Euclidean 'last' singular orbit 
becomes the sphere equator, which is also a straight line in the spherical geometry. Qualitatively, 
this reminds the results in the Kepler problem, where the single Euclidean parabolic orbit 'splits' 
for Ki < into a full family of Kepler orbits in H 2 , bordered by two different limiting curves, an 
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horoellipse as a limiting form of ellipses and a horohyperbola as a limiting form of hyperbolas, with 
a full interval of parabolas between them. The separating role in the oscillator problem is played 
by equidistants. 

Of course, since the nine CK manifolds S K1 r K2 i are geometrically very different, many dynamical 
properties display some differences according to some distinguishing properties of the manifolds; 
nevertheless, all these differences can be traced back ultimately to characteristics related to the 
signs of the basic parameters k±, K2 so there is a unique theory that is simultaneously valid for all 
the cases, for any value of the curvature and any signature type. 

The analysis of the orbits in the curved harmonic oscillator leads, in a natural way, to the 
theory of conies on spaces of constant curvature and any signature type. Although Sects. 3 and 4 
were mainly concerned with dynamical questions, Sect. 5 was written emphasizing its geometrical 
character. It is clear that the theory of conies on spaces of constant curvature is a geometrical 
matter of some importance deserving a deeper study within this CK formalism that we hope to 
present elsewhere; in the three Riemannian spaces of constant curvature conies have been discussed, 
at different depths and from differing viewpoints, in papers dating from a century or more (see for 
instance [36], p. 229 or [15]). Of course, none of these papers discussed the case of Lorentzian 
spaces, where it appears that this theory has been never presented systematically. 

Appendix: Geodesic coordinates on two-dimensional mani- 
folds 

Consider the generic CK space <S^j K2 ], for any values of «i,/C2. When K2 is positive it may be 
reduced to 1, and then this family includes the three constant curvature 2d Ricmannian spaces V 2 t . 
When K2 is negative, it may be reduced to —1 and the spaces are Lorentzian manifolds of constant 
curvature I? . We now describe the two types of coordinates employed in the paper. Choose any 
point O a point on S k j K2 \ and let l\ be an oriented geodesic (time-like if K2 is non-positive, generated 
by Pi) through O, and h the oriented geodesic orthogonal to l\ through O (hence space-like if K2 
is non-positive and generated by P2) (see Figs. 5 and 6) 

I: Geodesic polar coordinates 

For any point Q in some suitable domain (with time-like separation to O in the lorentzian case), 
there is a unique geodesic I joining Q with O. The (geodesic) polar coordinates (r, </>) of Q, relative 
to the origin O and the positive geodesic ray of h, are the distance r between Q and O measured 
along I, and the angle 4> between I and the positive ray l\, measured around O (Figure 5). These 
coordinates are defined in some domain not extending beyond the cut locus of O, arc singular at O; 
when K2 >0, <fi is discontinuous on the positive ray of l\, where there is a jump of 27r/y / Ki", but if 
K2 < then the range of <p covers all real line, and there arc no jumps, at the price of the coordinates 
being singular on the isotropes through O. In the sphere the domain of polar coordinates only fails 
to include two points: O and its antipodal; for the hyperbolic plane it covers all the space except 
the point O; for the Anti De Sitter this covers the domain with time-like separation to O, save O 
and its antipodal, etc. 

The expression for the differential element of distance dl is given by 

ds 2 = dr 2 + k 2 S 2 Kl (r) d<p 2 , 

so that in the standard Euclidean case (/ti=0, /t2 = l) we get cLs 2 | E2 = dr 2 + r 2 d(j) 2 . 
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II: Geodesic parallel coordinates 

For any point Q in some suitable domain, there is a unique geodesic V 2 through Q and orthogonal 
to l\, intersecting l\ at a point denoted Q\\ the geodesic V 2 is space- like in the Lorentzian case. 
Alternatively, for Q in some suitable domain, there is a unique geodesic l[ through Q and orthogonal 
to I2 intersecting I2 at a point denoted Q2', the geodesic l{ is time-like in the Lorentzian case. The 
points Qi, Q 2 can be considered as the 'orthogonal projections' of Q on the lines li, l 2 (recall that 
S Kl i K2 \ is in general not an affine space). Then we can characterize the point Q by (Figure 6) 

1. Two coordinates (u,y). The coordinate u is the canonical parameter of the clement in the 
one-dimensional subgroup of translations along l±, generated by Pi and with label Ki, which 
brings O to Q\\ this value coincides with the distance between O and Qi computed along 
li with the CK space metric. The coordinate y is the canonical parameter of the clement in 
the one-dimensional subgroup of translations along 1' 2 , generated by e uPl P 2 e~ uPl and with 
label kik 2 , which brings Q\ to Q; this value is related with the distance between Q\ and Q 
computed along l[ with the CK space metric by a factor y/1^2 (remark in the Lorentzian case, 
y is always real, yet both a space- like separation and y/li2 are pure imaginary). 

2. Two coordinates (x,v). The coordinate v is the canonical parameter of the element in the 
one-dimensional subgroup of translations along I2, generated by P2 and with label k±K2, which 
brings O to Q2', this value is related with the distance between O and Q2 computed along I2 
with the CK space metric by a factor ^/k2, so that v is always real even in the Lorentzian case. 
The coordinate x is the canonical parameter of the element in the one-dimensional subgroup 
of translations along l[, generated by e vP ' 1 P\e~ vP ' 1 and with label K\, which brings Q2 to Q; 
this value coincides with the distance between Q2 and Q computed along with the CK 
space metric. 

In the first case we have the parallel coordinates of Q relative to (O, l\) and in the second case 
relative to (O, I2) (Figure 6). In the (it, y) system the curves 'u = constant' are geodesies meeting 
ortoghonally the 'base' geodesic h, and the curves l y = constant' are equidistant lines to the base 
li, and intersect orthogonally u = constant. In the (x,v) system the curves 'w = constant' are 
geodesies and the lines l x = constant' are equidistant to Z2. Notice that in the general case, with 
non-zero curvature Ki ^ we have u^x and w^y; only in the case of flat spaces do the equalities 
x = u, v = y hold. 

The (u,y) and {x,v) expressions for the differential element of distance ds 2 are given by [37| 

ds 2 = C^ lK2 (y) du 2 + K 2 dy 2 , ds 2 — dx 2 + K2 C^(x) dv 2 , 

The three coordinate systems can be related by the general formulae of trigonometry in the CK 
space S 2 l[K2] [23]. 

T Kl (u) = T K1 (r) C K2 (4>) , S KlK2 (y) = S B1 (r) S K2 (0) , C Kl (u) C KlK2 (y) = C Kl (r) . 

In particular, the two parallel coordinate systems coincide when hl\ = reducing to ds 2 — 
dx 2 + K2dy 2 ; in the standard Euclidean case with K2 = 1 this reduces further to ds 2 = dx 2 + dy 2 . 
In this article we have made use of only the (it, y) coordinates, but all the CK formulation can be 
easily expressed also in the (x, v) coordinates. 
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Figures and figure captions 



V(r) 




Figure 1. Plot of the harmonic oscillator Potential V(r) as a function of r, for several values for the 
curvature. Upper curves correspond to positive curvature n = 2,1,0.5; the slightly thicker line to the 
Euclidean plane k = 0; the lower curves to negative curvatures k = —0.5,-1,-2. All the functions have 
the same quadratic behaviour around r = and the quadratic Euclidean function appears in this formalism 
as making a separation between two different behaviours (an infinite wall at finite r versus a finite plateau 
at infinite r). In the case K2 > the natural range of the radial coordinate includes only positive values, 
but when < positive and negative r appear naturally. 



Vef f (r) 




Figure 2. Plot of the effective potential V° (r) as a function of r in the standard Riemannian positive 
curvature case (ki = 1, k 2 = 1) depicted for several values for J . All these potentials are asymmetric wells 
with two infinte walls at r = and r = 7r/2. 
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Figure 3. Plot of the effective potential V c (r) as a function of r in the standard Riemannian negative 
curvature case («i = — 1, K2 = 1) depicted for several values for J . The angular momentum standard J x 
corresponds to the slightly thicker curve, where behaviour of the effective potential changes. Curves for 
values of J greater (resp. lower) than J x appear above (resp. below) this curve and correspond to an 
equivalent potential without (resp. with) a minimum. 



Figure 4 ABC. Harmonic oscillator orbits in a hyperbolic plane configuration space, depicted in the 
conformal Poincare disk model. Each figure displays orbits with a fixed value for the minor semiaxis b (or 
equivalently, fixed 'partial energy' E2) and several values for the major semiaxis a, ranging from a — b 
(circular orbit, in green), seven ellipses for increasing values of a (in blue), an equidistant curve for a — 00 
or a = 00 (in red), seven ultraellipses for decreasing values of a (in blue) and finally the straight orbit for 
a = (in magenta). From left to right, b is ranging from 'small' (Figure 4a) to 'large' (Figure 4c) values. 
The potential centre is at the origin, which is a centre of the conies. Colors have been chosen to represent 
particular and limiting conies: circle (green), equidistant (red) and straight line (magenta). For ellipses the 
pair of focus (not marked) are on the horizontal line; for the equidistant the foci are at infinity, as well the 
focal lines, which are orthogonal to the horizontal line h at infinity; for the ultraellipses one set of focal 
lines is orthogonal to the horizontal line. Notice only orbits with total energy smaller than Eoo intersect 
the horizontal line ii and come back to the initial point. Orbits with total energy larger than this value 
are not closed and go to spatial infinity. The two families in blue (ellipses and ultraellipses) are the two 
generic behaviours, as explained in the text. 
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Figure 5 AB. The 'polar' coordinates (r,<j>). The diagram depicts the geometrical meaning of polar 
coordinates (r, cj>) in a general CK space S%r i, both in the locally Riemannian case K2 > (left) and in 
the pseudo- Riemannian case k,2 <0 (right). In all cases h,h,l are geodesies, and Zi, h are orthogonal. The 
light cone through O is also shown in the lorentzian diagram. The coordinate r has label k± while cj> nas 
label K2- In the Riemannian case, the coordinate r is non-negative, only vanishes at point O, where polar 
coordinates are singular, and the angular coordinate (f> ranges in the interval [0, 2tt/^/~K2~] with the usual 
periodic conditions. In the pseudo-Riemannian case r vanishes along the isotropes through O (and would 
be pure imaginary in the shaded area with space-like separation to O); the angle 4> ranges in the interval 
[— oo, oo] and for a given <j> the natural range of r involves positive as well as negative values. 




Figure 6 AB. The 'parallel' coordinates (u, y) and (v,x). The diagram depicts the geometrical meaning 
of the coordinates (u, y) and (v,x), for the same situation and with the same conventions as in Fig. 5. The 
lines Zi,Za are geodesies through Q orthogonal to h,h respectively. The coordinates u,x have label m and 
are (locally) defined near O in both the Riemannian and pseudo-Riemannian cases. The coordinated v,y 
have label K1K2 and the corresponding geodesies are represented dashed; in the pseudo-Riemannian case 
this means these geodesies are space-like. In all cases the ordinary sign convention applies. When ki 7^ 0, 
1 / m and v 7^ y, and equality is a degenerate property of the flat case. See the text in the appendix 
for more details, and note that the natural interpretation of all coordinates is as canonical parameters of 
one-parameter subgroup of translations along the lines li,l2,l'i,l<2 or of rotations around the point O. 



